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1.1 SU(A/^) Heisenberg Models 

The use of large N approximations to treat strongly interacting quantum 
systems been very extensive in the last decade. The approach originated in 
elementary particles theory, but has found many applications in condensed 
matter physics. Initially, the large N expansion was developed for the Kondo 
and Anderson models of magnetic impurities in metals. Soon thereafter it was 
extended to the Kondo and Anderson lattice models for mixed valence and 
heavy fermions phenomena in rare earth compounds [1,2]. 

In these notes we shall formulate and apply the large N approach to the 
quantum Heisenberg model [3-6] . This method provides an additional avenue 
to the static and dynamical correlations of quantum magnets. The mean field 
theories derived below can describe both ordered and disordered phases, at 
zero and at finite temperatures, and they complement the semiclassical ap- 
proaches. 

Generally speaking, the parameter TV labels an internal SU(iV) symmetry 
at each lattice site (i.e., the number of "flavors" a Schwinger boson or a 
constrained fermion can have). In most cases, the large iV approximation has 
been applied to treat spin Hamiltonians, where the symmetry is SU(2), and 
N is therefore not a truly large parameter. Nevertheless, the 1/N expansion 
provides an easy method for obtaining simple mean field theories. These have 
been found to be surprisingly successful as well. 

The large N approach handles strong local interactions in terms of con- 
straints. It is not a perturbative expansion in the size of the interactions but 
rather a saddle point expansion which usually preserves the spin symmetry 

of the Hamiltonian. The Hamiltonians are written as a sum of terms Oj.O,,, 

''J ''J 

which are biquadratic in the Schwinger boson creation and annihilation opera- 
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tors, on each bond on the lattice. This sets up a natural mean field decoupling 
scheme using one complex Hubbard Stratonovich field per bond. 

At the mean field level, the constraints are enforced only on average. Their 
effects are systematically reintroduced by the higher-order corrections in 1/iV. 

It turns out that different large A'' generalizations are suitable for different 
Heisenberg models, depending on the sign of couplings, spin size, and lattice. 
Below, we describe two large N generalizations of the Heisenberg antiferro- 
magnet. 

1.2 Schwinger Representation of SU(A/^) 
Antiferromagnets 

The SU(2) algebra is defined by the famiHar relations [S",S^] — ieap-yS'^ . The 
spin operators commute on different sites, and admit a bosonic representation. 
Since the spectrum of a bosonic oscillator includes an infinite tower of states, 
a constraint is required in order to limit the local Hilbert space dimension to 
2iS'-|- 1. In the Holstein-Primakoff representation, one utilizes a single boson h, 
writing S+ = hW2S-h'<h, S- ^ {S+)\ and S"" = h^h-S, together with the 
non-holonomic constraint < h^h < 2S. The square roots prove inconvenient, 
and practically one must expand them as a power series in h^h/2S. This 
generates the so-called spin- wave expansion. 

Another representation, due to Schwinger, makes use of two bosons, a and 
b. We write 

S+ = a^ , S-=b''a , 5^ = i(a^a - , (1.1) 

along with the holonomic constraint, 

a^a + 6^6 = 2S', (1.2) 

where the boson occupation, 2S, is an integer which determines the represen- 
tation of SU(2). This scheme is depicted graphically in fig. 11.11 

There are three significant virtues of the Schwinger representation. The 
first is that there are no square roots to expand. The second is that the holo- 
nomic constraint (|1.2p can be elegantly treated using a Lagrange multiplier. 
The third is that it admits a straightforward and simple generalization to 
SU(A''). That generalization involves adding additional boson oscillators - N 
in all for SU(A^) - which we write as with /i = 1, . . . , iV. The generators of 
SU(iV) are then 

S,. - blb^ ■ (1.3) 
These satisfy the SU(iV) commutation relations 

The constraint is then 
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N 

E^K = ^6' (1-5) 

which specifies the representation of SU(iV). The corresponding Young tableau 
is one with n^^ boxes in a single row. 

1.2.1 Bipartite Antiferromagnet 

We consider the case of nearest neighbor SU(2) antiferromagnet, with inter- 
action strength J > 0, on a bipartite lattice with sublattices A and B. A 
bond (ij) is defined such that i €z A and j G B. The antiferromagnctic bond 
operator is defined as 

= '^^bj-b,a^■ (1-6) 

This is antisymmetric under interchange of the site indices i and j, and trans- 
forms as a singlet under a global SU(2) rotation. 

Consider now a rotation by tt about the y axis on sublattice B only, which 
sends 

a,^-b^ , b^^a^. (1.7) 

This is a canonical transformation which preserves the constraint (jl.Sp . The 
antiferromagnctic bond operator takes the form 



(1.8) 
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The SU(2) Heisenberg model is written in the form 
T~L = J ^ ^ Si • Sj 

= -^E(4A,-25^). (1.9) 

The extension to SU(A^) for iV > 2 is straightforward. With N species of 
bosons, (jl.Sp generalizes to 

N 

=EM.v (1-10) 

The nearest- neighbor SU(A^) antiferromagnetic Heisenberg model is then 

{ij) 

-^T.{T.srsT-Ns'\ (1.11) 

where ^ 

(1-12) 

are the generators of the conjugate representation on sublattice B. One should 
note that Ti. of p.lip is not invariant under uniform SU(A^) transformations 
U but only under staggered conjugate rotations U and on sublattices A 
and _B, respectively. 



1.2.2 Non-bipartite (Frustrated) Antiferromagnets 

For the group SU(2), one can always form a singlet from two sites in an 
identical spin-S representation. That is, the tensor product of two spin-S* 
states always contains a singlet: 

S®S = Q®l®---®2S . (1.13) 

For SU(A^) this is no longer the case. For example, for two SU(3) sites in the 
fundamental representation, one has 3 ® 3 = 3 © 6. One needs three con- 
stituents to make an SU(3) singlet, as with color singlets in QCD, and TV 
constituents in the case of SU(iV). This is why, in the case of the antiferro- 
magnet, one chooses the conjugate representation on the B sublattice - the 
product of a representation and its conjugate always contains a singlet. 

But what does one do if the lattice is not bipartite? This situation was 
addressed by Read and Sachdev [7] , who extended the Schwinger boson theory 
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to the group Sp(7V). This amounts to generahzing the hnk operator Aij in p.6p 
to include a flavor index m: 



N 



m— 1 
2N 

= E ^oa'^.a^,.' • (1-14) 

Here, the indices a and a' run from 1 to 2N. They may be written in composite 
form as Of ^ (m, /z), where m runs from 1 to and fi from 1 to 2 (or "f and J,). 
In this case, on each site one has — and = The matrix A^^, 
is then A^^^^^ = ^mn^^iv^ where e^^ — io^^^, is the rank two antisymmetric 
tensor. 

If we make a global transformation on the Schwinger bosons, with bia — *■ 
Vaa> bia' , then we find 

A, ^(c/Mc/)_,&,„&^.„, . (1.15) 

Thus, the link operators remain invariant under the class of complex transfor- 
mations which satisfy U^AU — A. This is the definition of the group Sp(2A^, C). 
If we further demand that [/ S U {2N) , which is necessary if the group oper- 
ations are to commute with the total occupancy constraint, we arrive at the 
group 

Sp(iV) = Sp(2iV, C) n U(2iV) . (1.16) 
For = 1 one has Sp(l) ~ SU(2). The particular representation is again 
specified by the local boson occupation, = ^\oPia- "^^^ Hamiltonian is 

Here, we have allowed for further neighbor couplings, which can be used to 
introduce frustration in the square lattice antiferromagnet, e.g. the J1 — J2 — J3 
model [7]. For each distinct coupling J^^ (assumed translationally invariant), 
a new Hubbard-Stratonovich decomposition is required. 

One can also retain the definition in (jl.lOp even for frustrated lattices. In 
this case, under a global transformation b^^ — > U^^ b^^, the link operator A^j 

transforms as Aij (U^'U)^,, b^^bj^, and invariance of A^j requires U*'U = 1. 
This symmetry is that of the complex orthogonal group 0(A^, C). Once again, 
we require U G SU(A^) so that the constraint equation remains invariant. We 
then arrive at the real orthogonal group 0(iV) = 0(iV, C)nSU(7V). For iV = 2, 
in terms of the original spin operators, we have 

-i4^..=-^i^.>l^.. (1-18) 
= - {S^S^ - SfS^ + S]S] + S^) {N = 2) . (1.19) 
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On a bipartite lattice, one can rotate by tt about the y-axis on the B sublattice, 
which recovers the isotropic Hciscnbcrg interaction • S^. On non-bipartite 
lattices, the N = 2 case does not correspond to any isotropic SU(2) model, 
and so one loses contact with the original problem. 



1.3 Mean Field Hamiltonian 

Within a functional integral approach, one introduces a single real field (r) 
on each site to enforce the occupancy constraint, and a complex Hubbard- 
Stratonovich field Q^jir) on each link to decouple the interaction. At the 
mean field level it is assumed that these fields are static. This results in the 
mean field Hamiltonian 

= ^ E \%\' + E (Qij 4 + Qtj Aj) (1-20) 

i<j i<j 

+E^^(^L^- - "J + (v^)"'/'E (-^-^^ + ' 

i ia 

where V, the volume, is the number of Bravais lattice sites, and where a runs 
from 1 to for the SU(7V) models (for which p = 1), and from 1 to 2N for 
the Sp(A^) models (for which p = 2). The field 0,^, which couples linearly to 
the Schwinger bosons, is conjugate to the condensate parameter which 
means 



d<l>*c, VVN 

Let us further assume that the mean field solution has the symmetry of 
the underlying lattice, and that the interactions are only between nearest 
neighbor sites on a Bravais lattice. Then, after Fourier transforming, we have 

^ ^ k,Q:,a' 

where z is the lattice coordination number, and where K^^, = 6^^, for SU(A^) 
and K^^, = A^^, for Sp(iV). We define 

4 = 7E'%«^'''' (1-23) 

^ 5 

where the sum is over all distinct nearest neighbor vectors in a unit cell. That 

is, —S is not included in the sum. The quantity 7]g = ±1 is a sign about 
which we shall have more to say presently. On the square lattice, for example. 
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— Vx^^''" + Vy^^'^^- For symmetric K^^^, owing to the sum on k, we can 
replace with its real part, while for antisymmetric K^f^ we must replace it 
with i times its imaginary part. We therefore define 

7k = - Y^'vs cos(k -6) if K = K'^ (1.24) 

^ s 

= - sin(k -S) if K= -K' (1.25) 

^ 5 



The sign ij^ is irrelevant on bipartite lattices, since it can be set to unity for 
all 6 simply by choosing an appropriate center for the Brillouin zone. But on 
frustrated lattices, the signs matter. 

It is now quite simple to integrate out the Schwinger bosons. After we do 
so, we make a Legendre transformation to replace the field </)j^ with the order 
parameter = {h^^)/\/VN, by writing 

G = F-Y, {^ic /3L + Aa) • (1-26) 
The final form of the free energy per site, per flavor, is 

BZ 

(1.27) 

where k = n^^/N, and E^^^ is the condensation energy, 

= A ^ |;3k, J2 + \zY, ^aa' {Q 7k /3k,./3-k,.' + Q* it /^k.a/^-k..') • 
k,a k,a,a' 

(1.28) 

The dispersion is given by 



+ -E^con > 



4>^^\zQlJ . (1.29) 



The fact that g is formally of order (assuming k is as well) allows one to 
generate a systematic expansion of the free energy in powers of 1/A'^. 



1.3.1 Mean Field Equations 

The mean field equations are obtained by extremizing the free energy G with 
respect to the parameters A, Q, and „• Thus, 

+ ' = ^/(S^ ("-^(^^ + ^) + E (1-30) 

BZ 
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BZ ^''^ 

= A/3k,„ + zQj^ Yl ^aa' P-Ka' ■ (1-32) 

a' 

Here, Tik(r) = (e'^t/-^ — l) ^ is the thermal Bose occupancy function. In 
deriving the second of the above mean field equations, we have also invoked 
the third. Assuming that the condensate occurs at a single wavevector k, 
the last equation requires that = at the ordering wavevector, ensuring 
gaplessncss of the excitation spectrum. When there is no condensate, /^^ „ = 
for all k and a. 

It is instructive to compute ( ) — {^r a^n a' ~ '^^aa' )' which serves 
as the local order parameter. After invoking the mean field equations, one 
finds 

{St ) = ^e^('''-'')-^/?L/3kw - E KcA'Kc' ■ (1-33) 

k,k' k,a" 

Note that the trace of the above expression vanishes on average {i.e. upon 
summing over R), and vanishes locally provided that the condensate satisfies 
the orthogonality condition 

E/5k«/?k'„=VE|/3kar- (1-34) 

a a 

In the case of an SU(A'') antiferromagnet on a (bipartite) hypercubic lattice, 
the condensate occurs only at the zone center k = and the zone corner 
k = TT. One then has 

( St ) = N{(31, /3^„, + /3o,,) e--^ . (1.35) 

Thus, Bose condensation of the Schwinger bosons is equivalent to long-ranged 
magnetic order. 

At T = 0, there is a critical value of k above which condensation occurs. To 
find this value, we invoke all three equations, but set the condensate fraction 
to zero. For the SU(iV) models, the minimum of the dispersion occurs at the 
zone center, k = 0. Setting Wj^^g ~ ^' ^® obtain the relation A = z\Q\. The 
first equation then yields 




(1.36) 



BZ 



For d = 1, there is no solution, and there is never a condensate. For d = 2, one 
finds Kj. = 0.19 on the square lattice [3]. Since k = 5 for the SU(2) case, this 
suggests that even the minimal S = j niodel is Ncel ordered on the square 
lattice, a result which is in agreement with quantum Monte Carlo studies. 
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Consider next the Sp(A^) model on the triangular lattice. We first must 
adopt a set of signs r]g. There are three bonds 2 3 P®'" ^^^^ along the 

directions aj^, a^, and SL-^ SL'2 5 where the primitive direct lattice vectors are 
a;^ = ax and aj = ^ax+^ay. Lattice symmetry suggests rj^ = +1, 773 = —1; 
and 7/3 = +1 (as opposed to all 77^5 = 1), resulting in [8] 

7k = i sin 6*1-1 sin 6*2 + 5 ^^^(^2 - ^1) > (1-37) 
where the wavevector is written as 

with 2 being the two primitive reciprocal lattice vectors for the triangular 
lattice. The maximum of |7i(.p, corresponding to the minimum of the disper- 
sion ojj^, occurs when 7,^ lies at one of the two inequivalent zone corners. In 
terms of the 9^^, these points lie at {61,62) = (^, where 7^ = — and 

at (6*1, 6*2) = (^, where 7^ = Sachdev [8] has found = 0.34 for 
the triangular structure. As one would guess, frustration increases the value 
of relative to that on the square lattice. On the Kagome lattice, which is 
even more highly frustrated, he finds = 0.53. 



1.4 The Mean Field Antiferromagnetic Ground State 



For a finite system (no long range order or Bose condensation) one can explic- 
itly write down the ground state of the SU(iV) Schwinger Boson Mean Field 
Theory It is simply the vacuum of all the Bogoliubov operators 



where, 
and 



/?ka = cosh6lk 6k„ - sinh6'k ^Lka' 



tanh 26k = — - 



A 



(1.39) 
(1.40) 
(1.41) 



The ground state wavcfunction can be explicitly written in terms of 
the original Schwinger bosons as 



Cexp 



|0) 



i5]e*^«tanh^k. 



(1.42) 
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For N = 2, using the unrotated operators and b\ the mean field Schwinger 
boson ground state ^E'"'' is 



,T,MF 



exp 



J2 



|0). 



(1.43) 



'^^^ contains many configurations with occupations different from 25* and 
is therefore not a pure spin state. As shown in [9], under Gutzwiller projec- 



tion it reduces to a valence bond state. Since tanhl 



k+TT 



tanh^. , where 



TT = (tt, TT, . . .), the bond parameters Uij only connect sublattice A to B. Fur- 
thermore, one can verify that for the nearest neighbor model above, u^^ > 0, 
and therefore the valence bond states obey Marshall's sign. 

Although 5'"'' is manifestly rotationally invariant, it may or may not ex- 
hibit long-ranged antiferromagnetic (Neel) order. This depends on the long- 
distance decay of u^j. As was shown [9, 10] the SBMFT ground state for the 
nearest neighbor model is disordered in one dimension, and can exhibit long- 
range order in two dimensions for physically relevant values of S. 

For further calculations, it is convenient to introduce the parametrizations: 



(20^^+2(1-7^) 



c i 
1 



^\Q\ 



T 



z\Q\ 



(1.44) 



Here, c, and t describe the spin wave velocity, correlation length, and the 
dimensionless temperature, respectively. In Fig. ll.2l the dispersion for the one- 
dimensional antiferromagnet is drawn. 

At the zone center and zone corner the mean field dispersion is that of free 
massive relativistic bosons. 



(20- 



ik-m2 



= 0, TT. 



(1.45) 



When the gap (or "mass" c/20 vanishes, cJj, are Goldstone modes which 
reduce to dispersions of antiferromagnetic spin waves. 



1.5 Staggered Magnetization in the Layered 
Antiferromagnet 

Consider now a layered antiferromagnet on a cubic lattice where the in-plane 
nearest neighbor coupling is J and the interlayer coupling is a J, with a ^ 1. 
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Fig. 1.2. Mean field dispersion uj^, in the domain — tt < fc < tt, for the one- 
dimensional antiferromagnet. 

We expect long range magnetic order at a finite Neel temperature Tn- The or- 
der parameter, which is the staggered magnetization M = ((—1)' e*'^'^ Sj^ j) 
becomes finite when the in-plane correlation length ^, which diverges expo- 
nentially at low T, produces an effective coupling between neighboring layers 
a^^(Tj^) which is of order unity. Here R locates the site within a plane, and 
I is the layer index. This means in effect that the coarse grained spins start 
to interact as if in an isotropic three dimensional cubic lattice which orders 
at Tn. The interlayer mean field theory, introduced by Scalapino, Imry and 
Pincus (SIP) [13] in the 1970's, can be appHed within the SBMFT. Here we 
follow Keimer et al. [14], and Ofer et al. [15], to compute the temperature 
dependent staggered magnetization, in the range T e [0,Tj^]. 
The Hamiltonian is given by 

^ = XI (^R., I ■ ^R+i, i + ^RJ ■ ^R+y, ; + " Sr.^ r Sj^ (1.46) 

The interplane coupling is decomposed using Hartree-Fock staggered magne- 
tization field: 

«5^,z5^,m^(-l)'e*'^-''(5R,;-^R,m) V/- — ' (1-47) 



Here h-^ ^ is the local Neel field due to any 
Assuming a uniform solution, ; = h, 



where it is assumed that M = M z. 

ordering in the neighboring layers, 
self-consistency is achieved when 

- = 2M{T,h) = 
a 



where M{T, h) is the staggered magnetization response to an ordering stag- 
gered field /i on a single layer. 
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T/T 

N 

Fig. 1.3. Numerical solution, from Ofer et al. [15], of the staggered magnetization 
M{T) of the layered antiferomagnet for various values of anisotropy parameter Ooff- 
Mo = S — 0.19660 is the zero temperature staggered magnetization, and Tjv is the 
Neel temperature. 



Extracting T-^ is relatively easy, since as T — s- T-^, /i — > 0, and the expres- 
sions for ( Oj^ ) and ( &r &r )■ In this limit, one finds that the second mean 
field equation in (|1.3ip is not affected by the staggered field, which simplifies 
the calculation. At Tj^ one finds 

2axlM = '^- (1-49) 

Since we know that xId ?2d(^n)' recover the ordering temperature of 
the SIP theory. The more precise calculation yields (restoring the Heisenberg 
exchange energy scale J), 

^ ^ _J_ ( 2^Mo 

N log a \\\og{\^nog{Aa/^)Kq)\ 

The numerically determined M{T) — h{T)/2a is shown in Fig ll.3l for various 
anisotropy parameters. 

One can also analyze the layered antiferromagnet using the SBMFT's na- 
tive decoupling scheme, without proceeding via the interlayer mean field the- 
ory of (|1.47p . Starting from an anisotropic Heisenberg model with in-plane 
exchange and interlayer exchange Jj_, one assumes a mean field solution 

where Q^j = Qy when (ij) is an in-plane bond, and = Qj^ when (ij) is an 
out-of-plane bond. The second mean field equation, (ll.31|l . then becomes two 
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equations. The Neel temperature can be written Tj^ = ,/|| f{Jj_/J^^), where 
f{a) is a dimensionless function. To find f{a), wc demand that the spectrum 
be gapless, but the condensate vanishes. This results in two coupled equations, 



J. J- 

^+1 = jd^\\ P||(7||)/d7^P^(7^) («(7||,7^) + i) (1-51) 

1 1 

/ dl^\ P\\ (7|| ) / d7± P±(7±) (7|| /f^) {n + \) 



9 T 1 1 

II ^i7,iP,i(7,i) /rf7^p^(7^)(7i/fi)(n+i) 



(1.52) 



where e = |(3_l/2(3||| and 0(7||,7_l) = ^^(1 + e)^ - (7|| + ejj_)'^ , and where 
= (e^/*<= — l) ^, with = T]vj/4 Once the above two equations are 



n 

solved for e and t^, we determine Qy from 
1 1 



g,, = J|| y d7|| P||(7||) / ^7^ pAi±) + ^) • (1-53) 



The functions Py and are given by 

^||(7||) = ;|k(1-7|) , P.(7j = . (1.54) 
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